NASA TECHNICAL MEMORANDUM 


NASA TM- 77381 

/U/9-S/? TM- 1738/ 

ir ‘ 


NASA-TM-77381 19840008422 


SHOCK CONDITIONS AND SHOCK WAVE STRUCTURES 
IN A STEADY FLOW IN A DISSIPATIVE FLUID 


Paul Germain and Jean-Pierre Guiraud 



Translation of "Conditions de Choc 'e^ Structure 
des ondes De Choc Dans un Ecoulement Stationnaire 
de Fluide Dissipatif", O f f i d£:cNa tional ; dA etude s et 
de Recherches Aeronautiques , Chatillon, -prance', Puhl. 
No, 1G5, IL962, 45 pages. HAM:->rc,->/v;na:MiA 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON D.C. 20546 DECEMBER 1983 


Utb 1 <J ' : T3 



NF00318 



STan&aHO titlk pack 


!• > »i rt N*. 

NASA TM - 77381 


4. Title «0»4 WVUtW 

shoc: 


7 » AwH 

Pau 


AmihIm H*.- 


X ■•dpi**!', C H». 


J. t n o *D 


t. P*.l 
SC 
Bo 


Tran 

Choc 

Offi 

Chati 


* 7 . K*j 


IT. Sh 

.Ua 


m 



,l- , oaoo 

■••'-j 

• ; 0".i:x> 

I'T'I 

• f * 

::o 



\r \ ©:ti \x\ xtsix-o. 
O iC :!>• © - :!: 2 l:tiO <!:u x.—- 

~s ■<: itci© or" m;:i> tnm 

I'd :d~ m i't'i to x- ■ ::u is: x- • — • m m 
o cr.i 'x: m ex i© i!D’ ©i to 1 o to •••- , • 

13 rj ‘ZCDCCIC: fr| -Cs. •• “j ■•:<!"CI 

r; OX rt'l 1-0 X'-ISiU ?x>haj 

iv i , :;o c;i c •— 1 3 : —i. ~i .-+• 


u. >; — !.li i:.i 

©to r.u :S to 

— | * O' — Mil I'd 

""i T^’lO . — i ■:'■ . — . 

to it' ■< 

o' < 

it' — I o 
ai in c o 
Q O Ip !< s " ijj 

O I'd <. —■ 

O'© tn n? cu : -- 
' ~"o •— c :p 

.-+■ ,-+0 .'P.'Tt ip 
:©•• — -’i; ,-i iX« 
it' O "s O X 
:J it' n 

o to O 111 111 
O " O Q. .“+• 

:p o 3 

< — lO O C It' 

if ©© c to 

ip.if ;;;i •* rs 


□Z:\ O'-H T| 
tne: -n — mu ■■■■' 

acor’ '©-.': 3 > 

;s::ot e>©' m © 

c:iti :;a :p 

33030 -N'.©'« 
COtO I'T'I 

n'i \ r:i to c” 

::?o —1:3:: 

X ■'■j'l 1 — 1 i o;i ■ 
Tiootnco 
com s: , ©r 

:x:co x i-~.i 

©to ;e:"ti — 1 

'• 0 C 03 >”i 30 
:n::o m. x- 
x iti 30 iti :p ;s: 
s <: 0 .1 o •* 
a:i' :x x 'f 

130 •—! si=: '• — O't 

3 Da>c:trj>.'iu 
. e.o o ■ x 

©I'—' I T 1 1 T j i— >.,.4. 

— 1 tP'tu 


T'| 

■ cox- 

D-. if 


!'j.i e 

:~i‘ roa, 
«»«•! • \ 


0 ~i cc« 

■ex;^ x. ■ / 

if I'J'l I — *•* / '—>■ 

• ;;k" •- ‘i:n 
•N.i.is.i-1 
O x.l'X 

■ O WOtO 

0 a't ^ "ti 

Q.. IT*' 

r: 

,'T- "'"ITI 

1 © o * — * c , *[ is:: 

id is: too x 

■ to —1 to v, ;s: 

4 i" c::ro ::D' 
111 " I'T'I X to 
:p' ' ::d- 


o 1 ti --i 

O Wf "vl 

I ?.' 1 I'J'l © to 

• X.s. I'T'I CCl 

:Ei IX'i 

Ill '-»■ 

< © 

If : 'S' 

1X1 1 -"' 


O 

1 — f" 

if 

c* 1 

© 0'’i 

Cl 

hi 

r.-0 


•• ,--l- 

c» 

::t:i 0 ■ 

<: ivuxi 


.H- 


O iZ. 



I'T'i :o 

i:u 


— 

i*— *• 

:r* io 0 

•—iif 

'O 

co::i> to ' 

T 1 -3 . 


" 'C! 

IVi ©1 

!U IQ •■1; 

;^r 

•”**5 ’ 

x col 

" .Is. i"i" 

•--* • 

O 

x ::d- 

— if Q 

if :;< 

O 


tom 

to 

«—+ 

O ■ — f 

tO if 


O 

— ■ :TI: 




m 

O , — 

13 * ■ 

m 

:x::s: :x> 

Ti m 


‘•■sj 

• !u 'l 

:p td 

O .H- 

— . 

1--' I'T'I ---I 


Ill 

■ if 

0 

if Ci.X 

tO if 


©1 1 


,-h . 

IXI 

1T1 70 

to 


O 

:r:::r: © 

O 



0 

*H-i . • 

13 

•~K 

••JfMl “SA 1 

is»., 

ai 

O 

...1 


,-t' if — 1*. 

-- • i;i;i 


I'T'll© 


’”?i 

w l 

•"Ot.O 

O C* 

:~j • 

h. 

iXiO 




X' ..lx 

:~3 

•a 0 

Q’ 1 

to:s: "t'i 



111 

\7".i 

1 V .c 

O -H- 

if 

v % s 1—1 r~ 

M 



I'T'I 

•<! 

"’’u 


O C.:« 


tu 

(Si 

CO 

m m ©■ 


:ti 

[ — 1 


1X1 

.-4 

. ::u 

tO if 

in 

if 

O' ’ V 



if 

"V 

•rz 

:p if 

•< 

-i:"D . 

— | 


m 

—I 

111 13 < 

.-h 


:x> m "T-i 

...... 

\Z' 

a 

c: •#! 

■-I- Hi Hi 

■H- 

O ■ 

! — ! — | — 

m 

IQ 


, u 

if '< rs 


■ 

:x> e: 


I“h. 


O 

— «. 

if "T. 

\Zi. 

I'T'I 1 1 — ■ 

CO 

0 

■H- 

j'”” irs!" 

,-y- O 

tO I'll 

110 

. 2 mO 





x- ::i> 

;:. 7 ' 0 e 

C .-f 


ITIO t .0 

.1 

•** 

O 

to ai 

If "T. tO 



pcizx 



•iid 

to X- 



jj-. 

!ki .1 

■ n*i 

c* 



'••T'I 

0 :d 

i * 1 * 1 ■ 

Cf 

::o © 

.-,•■ 4 * 


il!? 

1 — 1 ‘ 3 I| 

0:1 IXI 

to if 

if 

. to I'T'I ::t;i 

O' 

*• 


I'T'I 1 

to i-t- 

***r 





\;u 

O 

m 

m 


1 ' f‘h«— r if""* 

I'T'I 



•-X.I 

0 :p 

!H 

Iwf 

— ra i—; 




Iw» CO 

•H-.O 

rxi "r. 

111 

I'T'i ;ti ::&■ 

- icrt* 



0 CO 

if 

•< it' 

1:0 . 

30 tO 






CO 

If 

a;io n 

»■— 1 



if*""' 

©• ©' 


Cl, 

X X- '. 

tpi 



..iL 

if in 

if 


• — 1 

rj: 



iti :x 

. < - 

13 

O 

l—l 




:z: to 

if • 


*^i . " 

*w' 

■ 0 



— 1 


u 

MHN*. 


C«t»rW 


Tulas 


A I'is 4 , 306 





X OtnmMi AuiiiImH*.- 


1. R«mHN«, _ , 

NASA TM-77381 


4 . TtiU WUiI* . 

SHOCK CONDITIONS AND SHOCK WAVE STRUC- 
TURES IN A STEADY FLOW IN A DISSIPATIVE FLUID 


7 , 

Paul Germain and Jean-Pierre Guiraud 


ST AnOaRD TlTLK FAC ( 


X Recipient'* C*t«U( N*. 


December 19 83 


4. F*.f*rmlng Oi|mIi«Hm Cm4m 


I. Oi|*n<i«flM R*p*»« He. 


IS. WecV Unit Ns. 


t. futiminf Or|ani aetien Nan* m»4 A44r*t* 

SCITRAN 
Box 5456 


IX T ,p* *1 R*pad **4 Facial Ca*ac*l 

Translation 


12* Attncjr N*«i «i4 AWr*it 

National Aeronautics and Space Administration 
Waaningtoa, D.C. 20546 l' 4 * *H««r W. 


IX Saf planaatoy H»t*» 

Translation of "Conditions de Choc et Structure desondes.De 
Choc Dans un Ecoulement Stationnaire de Fluide Dissipatif , 
Office National d' etudes et de Recherches Aeronautiques , 
Chatillon, France, PuBl. No. 105, 1962, 45 pages. 


IX AlltfMt' 


This work differs from conventional studies on the same 
subject in that the thickness and camber terms, 
usually neglected, are here taken into account. More 
precisely, calling- e the reciprocal of the Reynolds 
number based on the shock wave curvature radius, the e 
terms of the first order are systematically taken ^to 
account. The most important result is a system of 
giving a correction of order e for the various RANKINE- 
HUGONIOT conditions. The suggested formulas may for 
instance have to be used instead of the conventional on 
to evaluate the loss of the total pressure across the 
tached shock wave which is found at the nose of a very 
small probe in supersonic flow. 


17. K*j W**4* (MnIN fry A«Ui»h4 


W. DllMWHMiWowl 


Unclassified and Unlimited 


IT. Wtrll, Cltlill, (*( A)* 
. Uacluaif 


X S*«K<*r Cliiiil. W Mi hi*I 

TV Mb. *1 Fa*** 

Unclaaalf lad 

60 


r\!'/s 4 , 306 














i 

SHOCK CONDITIONS AND SHOCK WAVE STRUCTURES IN A STEADY 
FLOW IN A DISSIPATIVE FLUID 


Paul Germain* and Jean-Pierre Guiraud** 


SUMMARY 


/ 3 ** 


This work differs from conventional studies on the 
same subject in that the thickness and camber terms, 
usually neglected, are here taken into account. More 
precisely, calling e the reciprocal of the Reynolds 
number based on the shock wave curvature radius, 
the e terms of the first order are systematically taken 
into account. The most important result is a system 
of formulas giving a correction of order e for the 
various RANKINE-HUGONIOT conditions. The suggested 
formulas may for instance have to be used instead of the 
conventional ones to evaluate the loss of the total 
pressure across the detached shock wave which is 
found at the nose of a very small probe in supersonic 
flow . 


INTRODUCTION 

The present paper will make a contribution to the classical 
shock wave theory. Strictly speaking, a shock wave in a perfect 
fluid is a surface along which the various parameters which 
characterize the flow (especially pressure) experience this con- 
tinuity. The shock conditions are relationships which relate 
the various discontinuities. Such a discontinuity surface in 
fact is only a perfect fluid model of a zone with a very small 
thickness, which exists in a real fluid which is slightly dissipative. 


-* r 

Professor at the Sorbonne , outside collaborator with ONERA 
**Research Engineer at the ONERA 
***Numbers in margin indicate foreign pagination 
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Witnin it these quantxties experience substantial variations but : 
these are continuous. The flow is always continuous, no matter 
how small the effects of the coefficients are which characterize 
the various dissipation mechanisms. In an exact sense, how can we 
legitimately speak of a shock wave in such a fluid? Is it possible 
to write down shock conditions which complement the classical 
conditions taking into account the dissipative effect? Such is 
the purpose of this paper. /4 

This has led to various studies which have now become 
classical [2, 6, 7, 9] and provide a valid model when the coefficients 
of the dissipation mechanisms are evanescent. Recently, the question 
has been discussed by Sedov-Mikhailova and Chernii [8] and also 
by Probstein and Kemp [5] , in order to improve the description when 
the dissipation mechanisms are weak and not evanescent. This was 
treated in a memoire which discusses a number of interesting 
questions. For reasons to be discussed below, the responses given 
by these various authors are not satisfactory. 

Our discussion will be divided into four parts. In the first 
part which continues in a sense our Introduction, we will analyze 
the work of Probstein and Kemp. In addition to its intrinsic 
interest, this will allow us to discuss the relationship between 
our problem and modern aeronautics. The formulation of the problem 
is then described exactly and intuitive analysis is performed in 
order to show the kind of results obtained and the terms which have 
to be implemented in order to obtain these results. The second part 
has the purpose of setting the simple properties of a class of 
functions, called N.S., which allows one to formulate in a mathe- 
matically explicit manner the term which is loosely called shock 
wave and shock conditions in a fluid which has dissipation mechanisms, 
for example, a fluid which satisfies the Navier-Stokes conditions. 

In the third part we will use a global method to formulate the shock 
conditions in the most general way. We will use a symbolic notation 
in the case of a steady flow. The calculations are made up to 
first order in •=— . The result shows that some of the additional 

KG 

terms which are introduced in the conditions (compared with the 
classical Rankine-Hugoniot conditions) require a precise knowledge 
of the development of the flow inside the zone in which the rapid 
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but continuous variations of the various characteristic quantities 
take place. 

The study of this "internal structure" of the shock wave is 
discussed in the last part. We indicate how one should set up 
differential equations which govern this structure and we proceed 
with the solution of these equations. Independent of the global 
method one can then obtain the shock conditions and one can prove 
the coherence of the general theory which is presented here. In 
this article we will limit ourselves to the consideration of 
first order terms, but verification of coherence was made for all 
orders. /5 

1. PRELIMINARY GENERAL REMARKS 

With the construction of artificial satellites and rockets 
having slow re-entry, there is increasing interest in the prediction 
of the aerodynamic characteristics of flying machines which fly 
at a very high speed and a very high altitude. Even though the 
stresses encountered and the heat flux received are small, because 
of the low atmospheric density, it may be important to take them 
into account for long flights. 

Two aerodynamic regimes can be considered to be well known 
at the present time, at least from the engineer-designer point of 
view. First of all we have the molecular regime and the continuous 
regime. The molecular regime is the regime for which the average 
mean free path is sufficiently large with respect to the characteristic 
dimensions of the obstacle, and for this the gas molecules can be 
considered as independent of one another. From the point of view 
of calculations of these forces the aerodynamic aspect of this 
regime does not introduce any difficulties. The only question 
which is not entirely clear is the mechanism of the reflection of 
molecules at the surface of the object. At the other extreme, 
the continuous regime is so well known that it is not necessary 
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to discuss it. 

Therefore, for the aerodynamic designer, the aerodynamic 
regime which is most important for extended flight at high 
altitudes may be the intermediate regime. This is a situation 
which exists when the flight is not high enough so that the 
rarefaction of the air is sufficient so that one can utilize the 
very siptple molecular regime properties, but is not low enough 
so that one could use these numerous results for the continuous 
regime . 

The intermediate regime which we will now discuss is 
characterized by the requirement for describing the microscopic 
structure of the gas medium and equilibrium. This is a structure 
which varies greatly from the Maxwell-Boltzmann structure, which is 
adapted to extreme regimes. The Boltzmann equation states the 
compromise which is realized between the tendency to move away 
from the Maxwell-Boltzmann distribution due to the motion of the 

molecules by inertia, and the trend to return to this cononical 
distribution due to the regularization effect of the collisions. 

In principle it should be able to predict the variation of the 
microscopic structure of the gas medium and consequently to determine 
the aerodynamic characteristics of the obstacles in the intermediate 
regime . 

In any case, the extreme complexity of the Boltzmann equation 
(integral-differential equation which is highly nonlinear) has 
discouraged all attempts to obtain explicit solutions which are 
different from the Maxwell-Boltzmann solution. Experimentation at 
high speeds and in a rarefied atmosphere is difficult. There are 
at least two decades of aerodynamics of the intermediate regime 
which are entirely unknown. 

* 

But the development in this area is rapid, and we now have 
both theoretical and experimental reasons which lead us to assume 
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that the variation of the aerodynamic characteristics is 
perfectly known from one of these extreme regimes to the other. 

If this is confirmed, together with a good knowledge of the /6 

continuous regime and the molecular regime, then approximation 
formulas will make it possible to predict the initial part of 
the variation from either one of the other two regimes. For the 
engineer, this will be sufficient information for interpolation 
with enough certainty of the variation of the aerodynamic 
characteristics in the intermediate regime. The various recently 
developed models may be sufficiently powerful in order to make 
a direct study of the intermediate regime and also to support 
interpolation methods. 

1.1. CRITICAL EXAMINATION OF THE MEMOIRE OF PROBSTEIN AND 

KEMP (P and K) 

Recently the study of heat transfer of the stagnation point 
was discussed in a memoire* of Probstein and Kemp [5] in order 
to extend slightly the results of the continuous regime already 
known into the intermediate regime. In reality, the work of these 
authors goes even further. Figures 6 and 7 of their memoire show 
this. It is possible to show that an interpolation within the 
intermediate regime between the continuous regime and the molecular 
regime is possible with minimum difficulty. Nevertheless, in spite 
of this technical success, and one might say even because of this, 
the work of P and K is not without criticism, both in the conceptual 
and fundamental sense. The purpose of the present article is to 
examine one of the points of the P and K memoire which can be 
criticized. We will attempt to shed light on the controversial 
question [4, 5, 8] of shock conditions within a dissipation medium. 

But before starting the present study, we have to somewhat analyze 
the work of P and K. After this is done, the reader will undoubtedly 
feel that the question formulated is a fundamental one. We are 
considering the problem of the flow over a sphere or a circular 
cylinder near the stagnation point. We will follow the investigation 

* In the following we will call this memoire by its initials P and K. 
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of P and K about the influence of rarefaction. One is then led to 
distinguishing less than eight flow regimes mentioned above. 

Regimes 1 and 8 are naturally the two extreme regimes mentioned 
before. The mean free path and the viscosity are directly related 
( n = p aA with y = viscosity and p = specific mass, a = speed of 
sound, \ = mean free path) . The continuous regime is a fluid regime 
which is almost perfect, that is, it has only slight dissipative 
mechanisms. It follows from this that the viscous fluid effects 
and the heat conduction are only felt in the boundary layer near 
the obstacle. Naturally they are felt through the shock waves 
because the concept of a shock wave, discontinuity surface are 
conditions of a perfect fluid, which have to be replaced by the 
path layer with a very high gradient if one takes into account 
the dissipation mechanisms. Also, the shock wave effects are one 
order higher than the boundary layer effects so that in regime (1), 
which one can call the limit boundary layer, the shock waves can 
be modeled by a discontinuity surface. The remaining six regimes 
establish a division which may be artificial of the intermediate 
regime. Using approximation formulas, one can distinguish in 
reality three sub-regimes within the intermediate regime. There 
is the almost continuous regime, or sliding regime for which one 
takes into account terms on the order of M„ //rTT (M„ : Mach 

number of the incident flow, Re M =p„u 00 R/fi„ : Reynolds number 

of the flow upstream based on the radius of the obstacle) . There 
is the almost molecular regime for which one takes into account terms 
on the order of Re~/M„ . There is the transition regime proper, 
for which one considers all of the characteristics of the intermediate 
regime. The regime (7) of P and K is the "First-order collision 
regime", but does not differ from the almost molecular regime. /7 

One also has to identify regime (2) or the "Vorticity interaction 
regime" and the almost continuous regime mentioned above. The re- 
maining four regimes: "viscous layer regime (3), "incipient merged 
layer regime" (4) , "fully merged layer regime" (5) , "transitional 
regime" (6) with the only real transition regime. From the point 
of view of kinetic theory of gases, this division into six sub- 
regimes is not necessary. 
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In fact, the point of view of P and K is different because 
it is not based on the concepts of kinetic theory of gases, but 
on the theory of continuous media and the Navier-Stokes equations. 

The best proof of this is that in regime (2) , in which we have found 
the sliding regime of the kinetic theory of gases, the sliding 
conditions and temperature jump conditions at the wall are not taken 
into account by P and K, and also not in regimes (3) and(4}. 

Since we have now discussed this, the rather artificial 
division of P and K is based on the substitution of the initial 
problem by four different problems which are assumed to be directly 
related to the evolution of the first problem as a functon of the 
degree of rarefaction of the air. Let us discuss this in more 

detail and also the four problems. The problem of regime (1) is 

so well known that it is not necessary to discuss. The problem of 
regime (2) is the same as for regime (l)but there is only a changed 
condition at the boundary layer edge, so as to consider the 
effect produced on it by an external vortex, a vortex produced by 
the shock wave. The problem of regime (3) is the problem of 

a flow of a viscous fluid between an obstacle and a shock wave, 

when one assumes that the Rankine-Hugoniot conditions are 
valid over the shock wave. The problem of regime (4) is the 
same as that for regime (3), but for different conditions over the 
shock, which makes it possible to take into account the thickness. 

Naturally the question can be posed as to how the four 
problems discussed are related to the evolution of the solution 
of the initial problem as a function of the degree of rarefaction 
of the air. This is a very difficult question and it does not seem 
possible to give a satisfactory answer, at least from the theoretical 
point of view. The only remark which we can make, and which was 
stated by P and K, is that recent studies have shown that the 
conclusions derived fronT the Navier-Stokes equations agree well 
with experimental results in the rarefaction domains, where from 
the beginning one would expect a good agreement. Based on this 
remark, we can replace the initial problem by a less difficult 
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problem and which is probably well posed: the variation of the 

flow of a viscous fluid over an obstacle having a rounded nose 
when the dissipation mechanism coefficients (viscosity and 
conductivity) increase from a value of zero. The question has 
special importance if one wishes to enlarge the area of debate 
and if one wants more than a simple analysis of an expedient means 
of predicting the aerodynamic characteristics at high speeds and 
high altitudes. In effect, this is one of the aspects of a 
fundamental problem which arises and which is being discussed more 
and more by applied mathematicians: how to study problems with 

nonlinear limits, in order to obtain explicit solutions and 
using asymptotic expansion. This is certainly the problem 
raised by the memoire of P and K, even though it is not strictly 
stated in these terms. It is necessary to expand the heat flux 
for example, according to powers of 1//Re„ , and the expansion 

starts with a term in 1 / /Re„ . The only valid criterion for 
distinguishing among regimes is the following: determination of 

each coefficient of expansion is related to solving a special 
problem, and one is certainly free to assign a flow regime to each 
one. Certainly, the preceding criterion is the only one that 
an applied mathematician can accept. 

1.2 EXACT DEFINITION OF THE PROBLEM BEING STUDIED /8 

Of course it is necessary to examine the memoire of P and K 
from this point of view. This is a problem of singular perturbations, 
because a unique formal development can only be valid close to an 
obstacle and also close to a shock wave. These are two regions 
where there are large gradients. Our purpose here is more limited 
because we will propose to only study the vicinity of the shock wave. 
But we will operate with a maximum of generality by posing the problem 
as a question in the fallowing form. When the fluid has dissipative 
mechanisms (essential viscosity and heat conduction) , which are 
small but not zero characterized by a small parameter e, and if we 
have a perfect fluid, then for this fluid, how can one develop 
using asymptotic expansions in e, a representation of the real fluid 
close to (E) , where (E) is the shock wave? 
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The answer to this question has to give conditions over the 
shock wave (Z) to be applied for regime (4). Conditions of this 
type were obtained by P and K but they are different from the ones 
established here and which are certainly exact, because they are 
obtained using an exact systematic procedure. The differences 
essentially come from the fact that P and K did not have a special 
technique for solving this problem. They wrote down the Navier- 
Stokes equations in integral form using intuitive ideas about the 
orders of magnitude which led to approximations. It is difficult to 
control the validity of these approximations and especially the 
significance of them with respect to the question being discussed. 
Again we have shown that obtaining terms of a higher order in an 
asymptotic expansion is most of all a technical procedure in con- 
trast to determining the first term. The first term often is obtained 
using intuitive methods. Only a posteriori will the result be 
justified by the technique. 

Let us now directly consider the question being discussed, 
by trying to develop as intuitively as possible the solution method. 

In Figure 1 a) we show ( 2 ) which is the discontinuity surface 
model of the perfect fluid. In a real fluid, it is part of a 
band B which lies between ( £ i) and ( 2 2 ) • Inside of this, the 
gradients are much larger than outside of it. In 1 b) we 
schematically show the profile of the densities when one moves 
perpendicular to (£) . Outside of B, we can make a formal expansion 
according to the powers of the small parameter e which characterizes 
the magnitude of the dissipation mechanisms. 



Figure 1. 
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Let 




p = p (o> + epf x) + .... 

+ € n l/ n> + 

p = p (o) + cp fl> + .... 

+ e“ p (o, + 

h = h <o) +ch (,) + .... 

+ € n h (o> + 

S = & o) + e^‘ ) + .... 

+ e° S (o> + 

v=V 0) +€V ( + .... 

+ c“V + 


( 1 ) /9 


be such expansions for the pressure, the specific mass, the 
specific enthalpy, the specific entropy and the velocity vector. 
When one limits oneself to variables with index 0, which is exact 
for a perfect fluid (e=0) and which constitutes the zero order 
approximation for the real fluid, the variables in question will 
satisfy the perfect fluid equation 

¥?+di v(p ( 0 V 0> ) =0. 

ot T 

w(°) 

p(0) m* g rad p <0)= 0' (2) 

0 

Dt 

and the Rankine-Hugoniot conditions over (Z) : 



[f/o> ^oj] 

= 0, 

[U (0 >] 

= 0 , 

[p (o, + p (o) w (o)Z ] 

= 0, 

[h«» + W< f ] 

= o . 


In the equation (3) the [ ] indicates a difference between 

the values of the expression inside the bracket, directly adjacent 

to either side of (E) ; w is the vector velocity component perpendicular 

* 

to (E) and U designates its projection on the plane tangent to (Z) . 
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We will designate the n order of approximation as being the 
result obtained by taking into account the first n terms in 
equation (1) . The approximation of rank n is formed by variables 
having an index n. The determination of the approximation of 

rank n includes the solution of a linear system of partial 

differential equations with suitable boundary conditions, 
especially on (E) . 

The expansions (i) are certainly not valid close to (E) , in 
the sense that they do not represent the real flow in this vicinity, 
that is, within the band B. Therefore, one has two expansion 

systems similar to (1) : one to the right of ( E 2 ) , and the other 

to the left of (£]_) , which respectively are their areas of 
representation. Even though the expansions discussed are not 
representative within the band, nevertheless they can be extended 
respectively to the left offJ^) and to the right of(Z^), and in 
particular, right up to (Z) , where they do not coincide. If we 
do not now discuss the question of which expansion is representative 
(1) close to (Z) , we can say that we have only a single system of 
expansion, valid to either side of (Z) , which experiences disconti- 
nuities through this surface. The question posed therefore is the 
following: How to express the discontinuities of the variables 

having index n through (E)? We will indicate such discontinuities 
in the following form 


[p (n) ]. [p (n) ], [h (n) ], [S (n) ] , [V (n) ] • (4) 

One would expect that knowledge of the shock conditions over 
(Z) ,that is knowledge of all of the discontinuities (4), allows one 
to determine each of the terms in the expansion one after the other 
(1) , using suitable boundary conditions at the various boundaries 
of the flow and the partial differential equation system which we 

mentioned. This then leads to a representation of the real fluid flow, 
which is valid except close to (Z) . The exceptional band becomes 
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wider as e becomes smaller. Inside this band it is necessary to 
carry out a special study. We state that one has to determine the 
internal structure of the shock wave. 

It is important to note that we are not especially interested 
in the vicinity of (Z); it is quite sufficient to have only the 
information provided by the expansions (1) . Then we have the situation 
represented in Figure 2 where (0) is the approximation curve of 
order 0, that is, the perfect fluid model, (1) is the approximation 
curve of order 1 corresponding to retaining two terms in the 
expansions (1) . (°°) is the curve obtained by retaining the complete 

series (1) , assuming that they converge. Therefore, it can be com- 
pletely legitimate to disregard completely the internal structure 
of the shock wave, but with the reservation that one must not lose 
sight of the fact that the fictitious fluid flow obtained is not 
representative for the real flow except for a more or less close 
vicinity of (Z) . 

However, it is necessary to 
determine the expansions (1) . For 
this purpose it is of fundamental 
importance to know the shock conditions 
for the approximations of various rank, 
that is, the expressions (4) for each 
value of n. For this purpose, we 
are led to pose the following question: 
is it possible to determine the discon- 
tinuities experienced by the character- 
istics of the fluid, for the rank and 
approximation, without taking into 
account the internal structure of the 
shock wave? We will see that the 
answer is negative and this is funda- 
mental because as a consequence, it is 
definitely not possible to disregard the 
internal structure, even if at the end 
of the calculation we are not interested in it, in contrast to what 
may have been stated above. In reality, for the perfect fluid model. 



Figure 2 
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that is for the zero approximation, it is not necessary to take it 
into account, but this is a very exceptional circumstance. Start- 
ing with approximation of rank 1, the determination of the shock 
conditions involves very greatly the internal structure of the shock 
wave. Sedov, Mikhailova, Chernii [8] nevertheless assumed the 
opposite and gave shock conditions with a first order of approxima- 
tion which were wrong. This error was corrected by P and K which 
were very well aware of the fact that it is necessary to take into /II 

account the internal structure of the shock wave using a fictitious 
expression for its thickness . But because there was no technique 
for the problem, these authors were not able to do this in a com- 
pletely correct manner . 

1.3 INTERNAL STRUCTURE AND SHOCK CONDITIONS 

Let us consider the one dimensional model shown in Figure 3 
which represents the profile of densities close to (E) ( z = 0) . 

The solid curve shows the behavior of a real fluid. The dash curve 
which is interrupted is the result of using expansions (1) . A 
study of the internal structure, therefore, amounts to studying 
how the solid curve can lead to a continuous agreement between the 
two parts of the interrupted curve in dashed lines . 

In order to study this internal structure, we are led to sub- 
stituting a variable £ for the variable z which is related to z 
with the relationship 

e Z = z, (5) 

if we note that the width of the band B in Figure 3 is effectively 
on the order of e . Let us now perform this variable change in the 
two functions shown in this figure: One corresponds to expansion 

(1) and is shown by dashed lines. The other one corresponds to the 
effective value taken on in the flow. Let us now assume that one 
could expand the result obtained according to powers of e and let 
us now consider in the two expansions terms of rank 0, rank 1... 
of rank n. The results obtained are shown in Figures 4 and 5 which 
we will now briefly discuss. Let us first discuss the nature of 
the curves obtained in Figure 4. When t, is finite, z is on the 
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order of e so that the approximation of order 0 based on (1) is 
independent of C to either side of (E) which results in Figure 4a. 

The approximation of rank 1 obviously is 
a linear function of £ to both sides of 
(E), and Figure 4b results. The approxi- 
mation of rank n is a polynomial of degree 
n in £ which results in Figure 4c. In 
this way, the curves derived from Figure 
4 are perfectly determined for any value 
of £ by the local behavior of the quantity 
under study and which is represented schem- 
atically in Figure 3 (p,for example) 
according to the definition given by (1) to either side of the value 
z = 0. Obviously, curves in Figure 5 are not very simple, but their 
relationship with the other curve corresponding to Figure 4 is 
easy to discern. When C is on the order of a few tens of units, 
the curves of Figures 4 and 5 must coincide for the corresponding 
cases. This is because outside of B, the deviation between the solid 
curve and the interrupted curve is exponentially small in e. Thus, 



Figure 3 
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for each order of approximation, the behavior close to z = + 0 of 
expansion (1) and the behavior close to K = + 00 of the internal 
structure are closely related^,- The approximation of rank n of the 
internal structure is equivalent to a polynomial of order n (up 
to an exponentially small term) , both for t = +cc and £ =-°° . The 
constant terms of these two polynomials are exactly the values 
taken on in (E) on either side by the approximation of rank n 
defined from (1) . We, therefore, have a rule which will be dis- 
cussed again, proven and used in the fourth part and which will 
provide a first method for obtaining the shock conditions . These 
will be a direct consequence of the internal shock structure. 

1.4 THE FIRST DEVELOPMENT OF A GLOBAL METHOD 


Disregarding any technical considerations, we would now like 
to discuss the principles of the global method which will be devel- 
oped in the third part. Without disregarding the intuitive consi- 
derations, at least for the simple case of the mass conservation 
equation, we will give almost the final formulation of the corres- 
ponding shock condition. As shown in Figure 6, we will consider a 
small cylinder perpendicular to (E) which cuts a surface element 
dA out of the surface which is limited by the closed curve (T) . 

By n we designate the unit vector perpendicular to (T) drawn in 
the plane tangent to (E) and which is directed toward the outside. 
The terminal bases of the cylinder with area dA are located at a 
distance of (E) which is sufficient so that the deviation between 
the approximation of order n based on expansions (1) and the beha- 
vior of the reel fluid are exponentially small in e. By w, we 
will denote the component of V along the cylinder generators. 

U = V-wZ is the projection on the plane perpendicular to the gen- 
erators. The distance to (E) along the generators (algebraic dis- 
tance) is naturally called z. Z designates the unit vector of 
one of these generators. The mass flux penetrating into the cylin- 
der through the two bases is - jpwdAj , where ( } designates the 
difference to be taken between these two bases. It is, therefore. 
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clear that the mass balance for the cylinder is written as follows 
where H is a geometric factor such that ds=Hds Q which is an arc 
element of the curve parallel to (F) to the distance z of (Z) : 


jp w dA | + £ ds„ J'^P U.n Hdz = o. (6) 

Now f designates an arbitrary characteristic of the flow. We 
call f its approximation of order 00 , which is a formal approxi- 
mation using expansions (1) . It is clear that according to the 
position of the terminal bases inside the flow, in (6) one can 
replace Jpw dA* by [pw dAj . After this, let us introduce the nota- 
tion ( ) + in order to indicate the difference between the terminal 
base and the section perpendicular on (Z) , from the side z > 0 or 
of the side z < 0. We then have 
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| p w dA | = [p w ] d A c + ( p w dA) + - (pw dA )_ • 


(7) 


After this, it is clear that even after the definition of the 
approximation of order infinity based on expansion (1) the 
following is true 

' ( p w dA)_ + f ds 0 j pu . n Hdz = 0 , 

r M (8) 

( (pw dA) + + / ds 0 / pU • H Hdz = 0 • 

T 


From (6), (7) and (8) and the remark made above that there 
is a time at which the shock condition on p w is the following, 
for order infinity 

[pw] = - 3 - 5— jr f ds °f (PU - pb) • n Hdz = o • 

The shock condition corresponding to the mass conservation 
was thus obtained up to r any desired approximation. Obtaining the 
other conditions does not involve any new ideas but the calcula- 
tion is slightly more complicated. Here we will not discuss this 
in this introduction, but we will give equation (9) a slightly 
different form. 
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Let us note that pu -pu is exponentially small in e when (z) 
satisfies both >> e and << 1. The points P and Q which intersect 
(£) at H in Figure 6 naturally satisfy this condition. But accord- 
ing to the definition of the approximation of order 00 for the inter- 
nal structure, we can modify the integral / by an exponentially 
small amount by substituting for P° its internal approximation 
of order % which is pu . But then £ = z/e is in fact the var- 
iable upon which the expression under the sign depends without 

these conditions, and except for an exponentially small error, the 
integral of (10) is equal to the following expression 

e (PU)* (M) = ey^"(pu - pu).n HdC- ( 10 ) 

Finally, the shock, condition (9) can be written in the form /14 


[pw ] 


lim 1 
d Ao - o d"A 0 


(p U )* • n ds . 


( 11 ) 


It should be noted that for the zero order of approximation, 
the second term is 0 and one finds as required the Rankine-Hugoniot 
conditions. It is clear that the second term will involve the inter- 
nal structure in an essential manner. When developing the internal 
structure as in (1) , we could develop (pu)* according to the powers 
of e and from this we could obtain the development of [pw]. It is 
then obvious that obtaining the shock conditions up to the approx- 
imation of rank n will involve internal structure up to order n-1, but 
not of order n. For the problem of P and K where one is restricted 
to the first order, it is sufficient to have the internal structure up to the 
approximation of order 0, which has been studied many times [2], [6], [7], [9] . 

It should be noted that the second term in (11) provides a local 
measure which is mathematically well defined of the shock wave thick- 
ness. There is a thickness for each conservation equation, that is, 
for each shock condition* In fact, (11) gives a displacement thick- 
ness because (10) is formally almost identical with the definition 
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thickness in the boundary layer theory. In this theory, the dis- 
placement thickness expresses the mass flux above the mass flux of 
the external flow which traverses a section perpendicular to the 
boundary layer. In the present case, the second term of (10) 
expresses the mass flux which traverses the lateral wall of the 
cylinder inside the shock wave, above the same mass flux relative 
to the outside flow. 

It remains to give an explicit expression for the second term 
of (12) by transforming the curvilinear integral into a surface 
integral. This is an easy calculation in vectorial analysis which 
will be done later on. 


The fundamental concepts and ideas which will be used during 
parts 3 and 4 have now been introduced in paragraphs 1.3 and 1.4. 
We will now give a precise definition of the functions which are 
involved in this work and will discuss a few properties which are 
particularly important for study. 

2. THE CLASS OF NS FUNCTIONS 
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2.1. DEFINITIONS 

Now we have to mathematically give an exact definition of the 
functions which represent the various characteristics of a viscous 
fluid flow near a surface which is a shock wave for the flow for 
the perfect fluid model. These functions depend on four variables 
x^, x 2 , z, e which have the following meaning: x 1 and x 2 are the 

curvilinear coordinates along the shock waves, z is the third coor- 
dinate which is 0 on the shock, z is a small parameter which char- 
acterizes the dissipation mechanisms. In this part, x^ and x 2 
will play the role of auxiliary parameters. In most cases, we will 
not use them for writing down these functions. 

The function f(z,e) depends on the parameters x^ and x 2 for a 
fixed z and small e. It behaves as a regular function of z and e. 
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The two functions f + (z,e) and f (z,e) obtained in this way for 
z > 0 and z < 0 obviously are not an extension of each other. See 
the discussion (1-1) . Also, in the shock zone which is studied by 
introducing the variable £ already defined in (1-5) we have 

e K = z. (1) 

Agreement between the upstream and downstream flow has to be imple- 
mented with this definition. According to already known results 
about the shock wave structure, this suggests defining the quantity 
as follows . 


DEFINITION 1 (class NS) 


The function f(z,e) may depend on parameters x^ and x 2 and 
appears in the class NS (feNS) if f is defined in the range 

K (a, e 0 ) : -a< z< a, 0 < e « e„ (a > o, e Q > 0) , and if in this range we 
can write f in the form 

f(z,e) = -||f + (z,e)+ f ‘(z,e)| + j -f(z,€)| tgh C+f (C.e) , (2) 

The functions f + and f appear in the class A, f in the class B 
which are defined above. We say that the functions (f + , f - , f) 
define f on K. 


DEFINITION 2 (class A) 


A function k (z,e) belongs to class A if it is defined and 
limited for any value of z when e«c 0 , and which can be developed 
in series of a power of e 


k (z.e) = 2 e ”k n (z), e < e 0 , (3) 

11 = 0 

The quantities k„ (z) are expandable in series of powers of z in 
the interval lz|< a , and f are zero outside a finite interval. These 
expansions are assumed to be differentiable an infinite number of 
times with respect to z, x^ and 
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Therefore, for the functions f + and f 
definition 1 we have 


introduced in the 
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( 4 ) 


f*( z ,e) =2eV (z) , 

n=o 

=2( ±Z ) P C . M <a - (5) 

P'0 

DEFINITION 3 (class B) 


A function f(C,e) is defined for any value of C when o<e<e„. 
For small e ? allows an asymptotic expansion 

f ~ £ c“4 (0 • (6) 

n=o 

For sufficiently large 1 1; | , f and f are exponentially 
decreasing (uniformly in e, as far as f is concerned) , that is, 
there are two positive constants A and a independent of e such 
that for example 

f(C,e)< Ae- a|<r| . (7) 

Finally, properties (6) and (7) are also valid for the deri- 
vatives of f with respect to x^, x^ ? . 

Two functions feNS which coincide in a domain K(a, e ) are 
considered to be identical. If (f , f , f) defines a function, 
(9/9/?) defines the same function f if f-g 

are zero over a domain K(a, e ) . It should be noted that (2) 
defines the extension of f for any ( 0 <e< € °) 

Now that we have made these definitions, we will give the 
exact definition of useful properties in the following. Here we 
will not again discuss the proofs in detail. 

2.2. PROPERTIES OF THE FUNCTIONS OF CLASS NS 

r 

We can immediately verify the following: 

LEMMA 1. If k(z,e) eA, g(C,e) eB , then 

k («C. e) % (C, e) . tgh C k (cC, e) g (C, e) 
are functions of class B. 
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and 


L 

.r t 


+ ^ 

From this we derive in particular that if f and fcA, f 2 (C,«) 
f 2 (C,e) eB, then 

—d* +n +-j(f - f) (tgh C + ? 2 K)) + f 2 (C) . (8) 

which defines a function of class NS. 

/17 

LEMMA 2. If f 2 (z, e) and f 2 (z,e) eNS, then there are also exists f^ + 
f 2> f^ f 2 , and derivatives with respect to x^ and x 2 of these 
functions . 


The only part of this statement which is not obvious is this 

state for the product of the functions.^ If we write 

f 2 - P 2 + Qj tgh C + f, , 

f 2 = P 2 + Q 2 tgh C + f 2 , 

f = f 2 f 2 = P + Qtgh C + f , 

then we have p _ p p + Q Q 


Q = Pj Q 2 + P 2 Qj . 

i = f 2 f 2 + f 2 (P 2 + Q 2 tgh C) + f 2 (P, + Qj tgh C) + Q, Q 2 (tgh* C - 1 ) . 
and the result appears as a direct consequence of lemma 1. 


LEMMA 3. (Outside development) 


For any z ^ 0 so that |z] < a, we have the following asymptotic 
expansions 


/ 


f (z, e)~ 




2 f n (Z) - if Z >0 , 


2 4' (z) , if z < o , 


(9) 


By definition, (9) defines the external expansion of f (z,e). 
We then have the hypotheses that the two series written in (9) are 
convergent within |z| < a and within them exactly define f + (z,e) 
and f (z, e) . 
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Then within K(a,Co) , we have the uniqueness of the func- 
tions t*, r, f which satisfy ( 2 ) for a given function f given of 
class NS. In effect, if f = 0, we find of necessity according to 
the preceding lemma that t* = f* = o and consequently f = o. 


LEMMA 4. (interior expansion) 


For any fixed ?, we have the following asymptotic expansion 

f (z,e)~ 2 *“ 4 W (10) 

n -o 

and the functions f are defined by 

n 

4 (O (p; (0 + 5 - (o)+y (p n + (o -p n - (o ) t g h c + ?„ (o , (11) 
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P n and are the polynomials in ? of degree n 

Pn 1 (0 =1 ( - C) P fi P , p . 

p.o 

The proof is immediate . 


( 12 ) 


LEMMA 5. 


For any function f (z,e) eNS, we can write down the identity 
f ( 6 ,e) -f(- 6 ,e) =V 6 P [f^ (e) +g(“. <0 . (13) 

p-O 

if o< 6 < a . In this identity g(C,e) is exponentially decreasing 
in > 0 ) and 


if 


[ f ]p ( £ ) = ^| e " (f".P * fn-P ) * 

n«o 

Proof is simple. In effect, we have 

f (6,e) - f (-6, e) = f*(6, e) - f'(-6,e) +g (-j, e). 

g (C.e) =-j - r «,<oj (tgh C - 1) + 

+ “]f*(-C,e) - r(-C,c)j (tgh? - 1) +f (C.O -f (-C.e)* 


(14) 
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The result is a a consequence of (4) and (5) and also of the 
reasoning given under lemma 1. 

LEMMA 6. 


If 6 is a positive number less than a, any function feNS can 

* 

be associated with a function f (e) such that 

fW)dz = c f* (e) + 2 £- P [F ]*,(«) +h 4.C) . 

p-i v * 

where h (C, c) decreases exponentially in C (C > 0) and 

[p Jp\i< e > = i cn ( c-i + f ; P -i > . 

n»o 

In order to prove this lemma, we can introduce a function 
f (z, e) eNS defined as the derivative with respect to z of 

\ Jf + (z.e) +F-(z,c) j +|- )f + - F- j tgh C , (17) 

This is a formula in which 


(15) 

(16) 



(18) 


Consequently, 


But we can write 


f*’ f (z, e) dz = 5 jf(z,e) - f (z, e)jdz +J^ f (z, e) 


r S f (z,e) dz = F + (6 J e) - F -(-6,e)+h 1 (±. e) , 
<5 £ 


dz 


^{f(z,e) - f(z,e)}dz =^Tf - ?) dz + h 2 (j,c) , 


hj (C.e) and h 2 (C,e) are two functions of C(C > 0) which decrease 
exponen tia 1 ly . 


The results (15) and (16) are derived directly 



( 19 ) 


ef*(c) =J + ” jf(z,e) -f(z,e)j dz * 
e/~jf (eC.e) - f(eC.e) [ dC • 

We can complete this result by stating that 
LEMMA 7. 


* 

The function f (e) which corresponds linearly with f(z,e) has 
the following asymptotic expansion for small e 

f* (e) ~ 2 e n f n * . (20) 

n»o 

with 

_p; (t> 

This result is obtained by writing formula (19) explicitly 

and taking into account the properties of functions of class A 

+ — - * 
for F and F and of B for f. We find that f can be calculated 

n 

from the n + 1 first terms of the internal expansion of f. 

In all the preceding definitions, the variables x^ and 
are introduced as parameters. In particular, [f] p (lemma 5) , 

[ F ] p , f*. (lemma 6) are functions of e and the parameters x^ and 

x 2 ’ 
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2.3. SHOCK CONDITIONS 

We can now indicate how the preceding concepts can be used in 
order to obtain the shock conditions which relate the functions of 
class NS which satisfy a relationship of the following form 


dt • (21) 


(1) Let us recall that f can be considered defined for any value 
of z when o < e < e 0 by (2). We can also verify that f*(f) only 
depends on values of f within K(a,e Q ) 
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( 22 ) 


(Af)** +^**bgdz = 0 • 

Here f and geNS, A and B are linear independent operators, 
independent of z which leave the class NS invariant. We find 
that the Navier-Stokes equations easily lead to relationships of 
type (22) . 


The results of lemmas 5 and 6 then lead one to write (0<6<a) 


[Af] o (e) + e(bg)* +g 6 P j[Af^ (e) +^[bg£ 1 (c)| = h (-|-,e) 


(23) 


The second term has an asymptotic expansion in e which is 
identically at zero for the asymptotic sequence of the e n . This 
is the same for the first one and consequently if we set (bg)* = 0 : 


[AfJ 0 + (b g) ;_ i+ V 6 p | [ Af n ] p + ~ [ bg n ] p :J=0. 


(24) 


The first term of (24) represents an expansion into a whole 
series of 6 which are identically zero. Therefore, all of the 
coefficients are zero. 


From this we derive that 

[Afl (e) +-[bg] p + . 1 (e)= 0 


and that 


y e n (bg)* =.V e" [ Af ] = - [ Af ] (e) = € (bg)* (e) 


We therefore can formulate the following theorem: 

THEOREM: If the functions f and g for class NS satisfy (22) for 

any pair of values z^ and then 


and 


[Afl (€) +e (bg)* (e) = 0 


[Af] p (e) +± I b g I*., (O s 0 


(25) 


(26) 


By definition (25) represents the shock conditions at z = 0 
which relate the quantities f and g which satisfy (22) . 
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h'e find that in this case, the asymptotic expansion of type 

* 

(20) of the quantity (bg) is a series expansion which is whole 
and convergent. 


The shock conditions of rank 0, l...n...are the conditions 
obtained by making the coefficients of power 0, l...n in e in 
equation (25) . 


[Af o ] 0 =0 

rank o 

iAfj o + (b g ); =0 

rank j 

I Af„ ] o + (bg^ j =0 

rank n 


3. DETERMINATION OF THE SHOCK CONDITIONS 
USING THE GLOBAL METHOD 


3.1 


>OTATION 


Che surface of the shock wave, a discontinuity surface of 
the perfect fluid model/ is called (I ) . On the surface we will 
consider a system of curvilinear coordinates x^, x^ for which the 
coordinate lines are the lines of curvature. By adding the dis- 
tance z to (E) measured along the normal perpendicular to it, we 
obtain a system of orthogonal curvilinear coordinates (x^,x 2 ,z) . 

We will assume that the lengths are measured with a certain unit L 
so that x^, x„ , z have no dimensions. In the same way, the pressure 
p, the specific mass p, the velocity V, the specific enthalpy h and 
the viscosity coefficient A,y and the conductivity coefficient k 
are assumed to be measured using the units p 0 = f>v o z , ^ ,v ,h o ,A o ,p o ,k o : 

We assume that p, p, h, y, y, k are dimensionless. The small 
parameter e which we have discussed will be 


E = 


Rc Po X, L 


( 1 ) 
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The unit vectors tangent to the coordinate lines will be called 
ei ,e 2 . e 3 and the elementary displacement will be called 

dM=hj e, dx, + h 2 e 2 dx 2 ' + e a dz , ^2) 

the square of the linear element is 

(3) 

da 1 = (hjdx ,) 2 + (hjdx 2 ) 2 + dz 2 • 

The partial derivatives with respect to x.^ and will be 
called 


_*L s f JL= f 

dx, 1,1 ' dx, ' 2 


(4) 


If we use N=e 3 for the unit vector normal to (£) directed 
towards positive z, we have 

N M = K 2 e, , N, 2 = K 2 e 2 , ( 5 ) 

where k 1 /H 1 ,k 2 /H 2 are the principal curvatures of (E) so that if 
an arbitrary point in space is written in the form 

M = m + zN - (6) 

the formulas (2) , (5) and (6) give 

hj = Hj + Kj z , h 2 = H 2 + K 2 z , (7) 

and which do not depend on x^ and , are the coefficients /23 

of the first quadratic form of (E) : 

dm 2 = (Hjdx 2 ) 2 + (H 2 dx 2 ) 2 • (8) 


The velocity vector y is written in the form 

V = u,ej + u 2 e 2 + we 3 • 


(9) 


3.2 NAVIER-STOKES EQUATIONS 


The unit tensor is designated by 1 and the viscous stress 
tensor is called et . The unit heat flux density is eq ; (S) 
is an arbitrary surface with an external normal N . Then the 
Navier-Stokes equations can be written in the form 
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(10) 





*« 

// 

“(S) 

IT 

*■ ( S ) 


p V . N dS = 0 , 

(pi + pVV + ct) . N dS = 0 , 
j ( h + — ) pV+e(t-V+q)|-NdS = 0. 


J 


In the case of stationary flow, this notation directly 
translates the laws of conservation of mass, momentum and energy. 


If D is the deformation tensor with a trace e = V.V 
we have 


t =-2 p D - A.G 1 

q=-k grad t, 

and one often sets 

, 2 

a = p v - y p . 


( 11 ) 


( 12 ) 


whare is the volume viscosity coefficient which is zero if one 
uses the Stokes hypothesis. 


Later on we will explain the calculations in the system of 
selected coordinates . But we will determine the shock conditions 
by using the dyad form shown above. 

3.3. SELECTION OF A CONTROL SURFACE. SYMBOLIC FORM OF THE 
SHOCK CONDITIONS 


The control surface (S) of equations (10) is selected in the 
fonn of a cylindrical surface as shown in Figure 7. The terminal 
bases are ( As ) and (A-s) . The lateral surface is formed by the 
generators of length 2 6 "perpendicular to (E) . They are supported 
at the centers by a closed curve (T) which delimits an area (A q ) 
°n (I) . 
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Figure 7 

Using the preceding notation, and when da designates the arc 
element of (F) , dS, the area element of (A + ^) or of (A q ) as the 
case may be, n is the normal to (F) , any of the formulas can be 
written in the form 

jf A • N dS - // A • N dS + / da n • f bdz = 0, (12) 

<A +S ) <A S > (D -S 

where A and b are either tensors or vectors. 

Equation (12) is of the form (2-22), so that the shock condi- 
tions are as follows in the symbolic notation and in integral form: 

// [A] -N dS + e f t> ■ n da = 0 . (13) 

J u o ) 1 in . 

In order to obtain the local shock conditions from (13), it is 
sufficient to transform the curvilinear integral into a surface 
integral, so that the first term in (13) can be written as a single 
integral taken over (A q ) . Application of a classical analysis 
theorem allows one then to conclude, since (A q ) is arbitrary, that 
the quantity to be integrated is identically 0 and it is conti- 
nuous in x^, x^ according to hypotheses. In order to perform this 
transformation, it is useful to distinguish the case where b is a 
vector or where b is a tensor. 
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First case: b is a vector: 


We can set 


and also 


therefore 


B = b* e 3 + b* e 2 + b* N , 

A = a 2 e x + a 2 e 2 + a 3 N , 

n da = Hj e 3 dx 2 - Hj e 2 dXj , 
f b-nda = f b* H 2 dx 2 - b* H t dx 3 


(D 


(D 


■•LKn- 1 .. I hT^ 


( 14 ) 


(15) 


For these conditions, (13) leads to the shock equation 

= 0*. (16) 


[A >] + h^T 2 | 'b.HJ./O.H.), 
Second case: b is a tensor: 
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Using the summation convention for silent indices, we can 
here set , 

b = btj ei ® ej . A = a„ e ; ® e, ; (17) 

Formula (15) must therefore be replaced by the following: 


/ B • n da = jj‘ j(b* 1 H ej) + (b* H e ) ( j—- (18) 

<n (a 0 > ' - 1 » ) n i n 2 

Let us recall that with the present notation, the expressions 
for the first derivatives of the vectors are as follows: 

b 12 b, 2 

e t 4 = — i- e 2 - K t e 3 . e = e, , e = K 2 e 2 . 

V h 2 

b 2 j bj j 

e I2 - — e 2 . e 22 e, - K 2 e 3 . e 3 2 = K 2 e 2 . 

hj hj 


6i j “0 * e 2 3 ~ 0 , 6 j,j 


= 0 
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Under these conditions, (13) leads to the scalar shock equa- 


tions 


/ e 

[A u ] + + (bj 2 Hj)^ ^ + b* 21 H 1>2 + b* H 2 K, - b * 2 H^j 


H.H, 


h,h 2 


[A lJ ]+; + Hj) -b;,H + b* a H 2 K + H, , 

w u ( » 1 ** ^ ' 


0 . 


= 0 . 


(19) 


[A 1S ] + I (b*j Hj) + (b^HJ - bJjHjK,- b t * 2 H 2 K 2 

h 2 h 2 ( * '* 1 
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3.4. SHOCK CONDITIONS 


It is only necessary to use the preceding results (16) and 
(19) in the Navier-Stokes equations written down as (10) . We will 
now examine the case of the mass equation, the momentum equation 
and finally the energy equation. 

3.4.1. Mass equation 

This case is particularly simple because A and b are 
vectors , , 

h 2 h i . 

A = p V , b=puj — e 2 + pu z -jj- e 2 , 

and (16) shows that, considering (9) , that 

( (pu,-^)* H 2 ) + ((pu.f^l’H,) • ( 20 ) 

This is a general equation. If we set e = 0, this gives the 
zero order condition [pw] = 0 which relates the values of p w in 
the external approximation of order zero. This is one of the 
usual shock conditions of Rankine-Hugoniot in a perfect fluid. 

9 

If as we have done here, one retains the terms of order e, the 
first term has to be evaluated with the two first terms of the 
external expansion of p w. The expression in ! then only in- 
volves the first term of the expansion of each of the starred 
quantities, that is, the term which is independent of e in the 


[pw] 


HjH 2 
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expansion mentioned. According to lemma 7 of part 2, this term 
can be determined from the first term of the internal expansion 
of the corresponding quantity. The classical theory [2, 6, 7, 9] 
of the shock wave structure has exactly the intent of calculating 
the first term of the internal expansion. In particular, we can 
see that (u^^ and (u 2 ) q are independent of X and respectively 
equal to U 1 and U 2 , which are tangential velocity components of 
the flow evaluated on (£) in a perfect fluid. Consequently, the 
shock condition translates the mass conservation as follows for 
the degree of approximation considered: 

['"']■ -hx I + • (21) 


In order to remove any possible ambiguity, let us discuss 
this equation in detail. The comments given for this simple case 
are also valid for other shock conditions which will be obtained 
in the following two paragraphs. The external expansions of p 
and w are 

P +ep + , z > 0 

p (o) - + c p (l >- + z < 0 (22) 



+ 

o 

1 

+ 

c w (1)+ + 

i w (o) - 

+ 

c w (1 >- + 


z > 0 
z < 0 


from which we find the external expansion of p w: 



( 

p(o)+ w <°> + + € 

(p<°> + 

w UH + pU> + 

w (o)> ) 

+ .... . Z 

> 0 

p w 

= ) 
l 

p (o) - w (o) ' + e 

(P (o) - 

w ( ‘V +p <l>- 

w (o) - ) 

+ .... , z 

< 0 

have 








[p w ] 

= P ( 

» )+ w (o)+ - p (o) " 

w (oV 

+ 





+ 

c j p <°> + 

-p<°>- 

w (1) - +p (I)+ 

w (o,+ - 

pUV w (o>- | + 




The various expressions which appear in this relationship are 
functions of x^ and x 2 . In this way, we have explicitly written 
down the first term of (21) . We will now explicitly determine 



the second term, that is essentially 
second part, we have 


p(x 1# x 2 ) = p (o) (x,,x 2 ) +0 (€) 


* 

P ^ 



According to the 


= / + ~ j P ' (o) ( C ) - j (P (o,+ + P (oV ) - \ (p <0)+ 


.^^(tghc+g^) ;dc+o(c), {23) 


where p° (C) is the first term of the internal expansion of p 
defined by 

p (x t ,x 2 ,eC.e) = P (x 2 ,x 2 ,C.€) = p (o) (x,, x 2 , C) + 0 (e) • 


The function p (o) (Q is precisely the one given by the classical 
studies on the internal structure [2, 6, 7, 9]. Since it satisfies 
the condition 

t (o ) . a. . _ 

P (*<=)= P 

the integral of (23) certainly is convergent. 


With this precise discussion, formula (21) is written 
explicitly in the following form 

I p(o>+ w (o)+ _ p Co)- w (o)- _ 0 

HjH 2 S p (o)+ w (n+ -p (o )- w (1) ‘+ p (1) V o)+ -p <1) * w ( o) - ( = 

1 » (24) 

= -^7 (p* <o) U jH 2 ) -- 5 |_(p^>U 2 H 1 ) , 

and the various expressions in this formula are obviously functions 
of x^ and but not of z and of e. 

3.4.2. Momentum equation 

It is now convenient to apply (19) with 

A= pi + pVV + et , b = A ej ® e t + -^ e 2 ® e z ) ; (25) 

but before giving the explicit form of these equations, one must 
determine how t depends on the derivatives of the velocity compo- 
nents. In this system of coordinates which we are using, we have 
according to (11) 
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( 26 ) 


with 


t = tjj ei ® e; . 

t ;i *> - 2y Djj - * fl6 ii , 


u i,i u * h >,2 


'11 


u. 


w 

u. h. 


K,w 


K 2 w 


*2,2 -*l l *2,l 

D = + + 

" h 2 h,h 2 h 


b w 


b z 


2D 21 = 2 D j2 =- 


2. 1 


* 1,2 


U l h l,2 

hjh 2 


U 2 h 2,l 

h l h 2 


(27) 


2D, 


2D 32 = 


6u 2 w j2 K 2 u 2 

6z h, h 


i)u 2 w 


i ^i u i 


2D 31 = 2 Dj 3 = 


bz h. 


6 = D., + D„ + D, 


ii 


If we only consider terms of order e in the shock conditions 
as we will do here, then in order to calculate the term in e in 
[ £ t] , it is appropriate to only consider in (26) and (27) the 
terms which are independent of e, that is to calculate (26) and 
(27) with the values of u^, u 2 , w (and possibly T if we assume 
that X and y depend on temperature) in a perfect fluid flow. With 
this rule we can explicitly write down the [A.,,] in (19) . It is 
still necessary to find r (b^) explicitly and it is sufficient to 
do this with the approximation of order zero. There are three 
categories of terms: p*6ij,(puiUj)* , (et;,)* . The terms of the first 

category can be immediately written down. Among the terms of the 
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second category, one can distinguish three types: 

(pu*)* = puj , (pu lUj )* = ^U t U t . Finally, ( P wu t )* = U t (pw)*=0 (these eva- 

luations are based on the fact that u^, u 2 , p w are constant in 
the internal structure of order zero) . Among the terms of the 
third category, it is only necessary to retain those of order zero. 
For this purpose, one writes down the internal expression etij 
of €tij and one only retains terms which are independent of e. 

We then find: 


^ + 0 (e) = et zz , 

eT„ = - (X+2P)-|j + 0(e) , 

= ‘ ** "5^ + 0 (e > = 0 = €T 32 » 

eT » = - ¥ ^r + 0 < c > = 0 < e ) = • 

et iz = 0 (e) • 


(28) 
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Therefore, one only has to retain 


“(3 |f)’ 

((«„)''” -((X+2ii| 


In these expressions -gp designates the derivative of w 
expressed as a function of x,, that is, the derivative of the func- 
tion which gives the velocity profile in the classical internal 
structure of the shock wave [2, 6, 7, 9]. 


Considering this, the shock conditions are easily written 
* 

as follows : 


* 


There is a sign error in formula (7) of reference [1] . This 
error is corrected in the present formula (30) . 
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I 


Ip 


wu t ] = e 



1 dw 

+ — (X p)* - 

K ' l 


1 } q \ 

— (p ‘ u * lH *>., + (P U iU 2 H,) 2 +p *u 1 u 2 h 1>2 - pu*H 2<1 , 

n i n 2 L) J) 

a) 



W » 2 
L 


1 1 ^ dz 

r 

H, 

H, il 


Hj d C 




h»h 2 


[p + pw 2 ] = e 


^ IP + IP 

r. ‘h 1 1 + 

J ' W W ' tin JJ 


(30) 


dw 

dz 


\ 


( A+2p) 4- X w ^ H I + 

dw \* / K, K 


H, H 2 


■K.u! K,U 2 


+ P (— ♦ _) t (p.x^)-(± + fi) . 

v h, h. / v d< ; Vh, h 2 M 


c) 


Let us repeat that the first terms have to be calculated 
with the external expansions for the approximation of order 1. 

The second terms have to be calculated with the internal expan- 
sions at the approximation of order zero if these are terms with 
asterisks. This has to be done with external expansions of the 
approximation of the order zero also for the terms [ ] . In order 
to remove any ambiguity, let us state that 



(31) 
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which corresponds to the general model (23) considering the fact 
that 



3.4.3. Energy equation 


36 




This time one has to apply (16) with 

A= (h+\) pV +e (t.V+ q } . b= A.(|i e , ® e i + ^- e , ® e *) . (32) 

but since the mechanism is now well known, we can proceed very 
rapidly. Let us only say that the term £ q gives a contribution 
of order e to the discontinuity terms (term between [ ]) . This 
contribution is calculated using the perfect fluid values. Con- 
versely, there is no contribution to be retained in (eq)* because 

* 

e£ji = °( € ) and e ^ 2 = °( e ) according to (11) and because (b^) does 
not occur in (16) . Without any difficulty, one obtains the fol- 
lowing result: 
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3.5. CONCLUSION 
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The formulas (21 ) , (30) and (33) translate the shock condi- 

tions when one takes into the account the terms of order e and 
only those, and the higher order terms are ignored. Let us recall 
their interpretation. If e is assumed to be zero, one finds the 
usual Rankine-Hugoniot conditions. If e is not zero, these for- 
mulas, therefore, give a correction to the R-H conditions, a 
correction which is of order e = g- . The Reynolds number is cal- 
culated with a length related to the average radius of curvature 
of the shock wave. The first terms contain the discontinuities 
of certain expressions which have to be evaluated with the terms 
of the external expansions up to order e inclusive. In the 
second terms, we have, with a factor of e, on the one hand the 
discontinuities (terms in [ ]) which have to be evaluated up to 
the approximation of order zero for the external flow (liquid 
fluid) , and on the other hand, the terms with asterisks which are 
calculated from the results of the study of the classical internal 
structure [ ] (first term of order zero of the internal expansions) . 
More precisely, there are four such terms with asterisks: 



(34) 


which have to be evaluated as a function of x^ and x as integrals 
in C over the interval - 00 , + 03 , of functions of x^, x^, £ . The 
general calculation formula is 


f*(Xl> X 2 ) / jf (0) (3C,.X 2 ,{;) (f M *(X lt X 2 ) + f <0> '(x i ,X I ) ^ 



(°)+ 

f (Xj, 


*.) 


- f (o) (Xj,x 2 ) ) 



dc. 


(35) 


and we find that for 



the terms 



become zero. . 
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The formulation (21) , (30) and (33) of the shock conditions 
is essentially related to the selected system of coordinates. The 
symbolic formulation (13) is intrinsic on the other hand. If it 
is advantageous for other reasons, there is no disadvantage for 
performing calculations in another system of coordinates on (E). 

For this purpose, one will write (13) by taking an elementary dS 
for (A q ) , for example, limited by an infinitesimally small parallel- 
ogram which is constructed using coordinate lines on the surface. 

The term [ A].N dS is immediately interpreted as the variation of 
the flux of a mechanical quantity (mass flux, momentum 

flux, energy flux) through the area element dS . This is the var- 
iation which is experienced when one passes through the shock waves, 
more precisely the discontinuity surface (E) . The term (A 0 ) 
obviously translates what escapes along the edges of the 

area (A ) within the thickness of the shock wave, beyond what 
o 

escapes if one were to use the values of the perfect fluid model 

to either side of it. This excess is precisely translated by the 
* 

operation . In the calculation of this term we have all of the 
terms involved in the Navier-Stokes equation and the curvilinear 
integral expresses an average value with respect to z of the 
various contributions of the flux through a cylinder with a base 
(F) and with a height e perpendicular to (E) . From this we obtain 
the following intrinsic rule: 

e / b* . n 

Id 

The variation of the flux density of a mechanical quantity 
(mass, momentum, energy) when passing through the shock wave (E) 
is equal to the sum of the two contributions if one ignores the 
terms of order equal to or greater than 2 in e . The first contri- 
bution is the one from viscous stresses and the heat flux which is 
adjacent to the shock wave on either side, and which is evaluated 
in the immediate vicinity to either side of the shock wave with 
the data supplied by the perfect fluid model relative to surface 
elements along the two faces of (E). This first contribution auto- /32 
matically is of order e. The second contribution can be evaluated 
formally by proceeding just like when one establishes the Navier- 
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Stokes equations from the integral form using the method of 
elementary parallelepiped volume- elements limited by coordinate 
surfaces. This contribution corresponds to the lateral part of 
the elementary volume surface. This lateral portion is a cylinder 
perpendicular to (I) and has the height e. However, the follow- 
ing precautions have to be observed. 

1. In the stress terms and the heat flux terms, only those 
have to be retained which are associated with a gradient perpen- 
dicular to (£) by dividing them by e; 

2 . replace the integration along the generator of the cylinder 

by an integration of an average value obtained by first performing 
★ 

the operation . In other words, one uses the product of this 
average value and e, and the contribution obtained in this way is 
of order e. In this calculation, it should not be forgotten 
that the operation * has no effect on the tangential velocity com- 
ponents. This rule is the natural extension of the corresponding 
rule which is applied when one establishes the Rankine-Hugoniot 
conditions using the method which can be "physical" and which is 
used in all elementary courses. This is a rule which was used by 
Sedov, Mikhailova and Chernii [8], but only the first contribution 
is shown there which gives an erroneous result. Probstein and 
Kemp [5] were not successful in giving a satisfactory form to the 
second contribution. 


r 
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4. STUDY OF THE INTERNAL STRUCTURE 


/33 


The results given before show that for determining the shock 
conditions of order 1, it is necessary to know the internal struc- 
ture of order zero and also the shock conditions of order n 
asstime that the structure is known up to order n - 1, even though 
explicit calculations have not been made. In this part, we will 
show how one can form differential systems which allow one to 
study the various approximations. In order to simplify the anal- 
ysis, we will only give the explicit calculation in the case of 
the order 1 structure. We will discuss the main properties of 
the differential system formed in this way and we will show that 
the internal flow "agrees" well with the external flow. One of the 
consequences of this study is to provide a new method, independent 
of the one which was discussed in the third part, and which gives 
the shock conditions of order 1. 

4.1. EQUATIONS OF THE INTERNAL STRUCTURE OF ORDER 0 AND 1 

The Navier-Stokes equations can be written in the form 

-* 

d i v (pV ) = o , (1) 

grad p + div (p V V + et)= o , 

div ( jh+^j pV +C jtV+q' ) =0 • 


This involves the derivation operator (nabla) which is given 
by the following using internal variables: 


± e 3.i_ + 

e 6C 


_ 6 _ 

dx. 


e t a + e. 


a 

dx. 


(2) 


so that symbolically, each equation (1) can be written in the form: 


= (e 
ac U 


a A + e^ 

dx, h. 


aA \ 

ax 2 / 


(3) 
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A designates vector or a tensor of order 2 , A is the same 
vector or teasr* with internal variables . 

If in Aone separates the terms which contain e as a factor 
from those whh do not contain, by writing 

A = B + cl . (4) 

and if we nofethat e, , e 2 » e are independent of C we can inte- 
grate this egfcion and we find that 

/34 

e 3 • B (C) = ~ e£ (A) , ( 5 ) 

with 


•^A 

a A 

6x i 

k + 


-A 

aA 

dX 2 

K • 

(6) 


Equation'S) is the symbolic equation of the internal struc 
ture. From dSs equation, we formally derive the equations of 
the internalatfcucture of the various ranks by setting 


A = 

oo 

2 

n = o 

-w 

c" A , 

B = 

oo 

V 

J-i 

n = o 

-(n) 

e" B , 

I = 

oo 

2 

_ (o) 

e n L 

£ = 

oo 

2 

-(d) 

e n £ , 


and by equatar the coefficients of the various powers of e we 
find 

e, • r ■ - jT' ; (/•’.<>). (7) 

We mustaas give the explicit form of the calculations. We 
will only dorBis for the terms of rank 0 and 1. 

4.1.1 5ks equation 

In this33tse A and the vector pV and L is zero. Con- 
sequently, (5 is written as 
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P w + £ Co * 0 , 


( 8 ) 


C Q designates a primitive with respect to t, 



as one can see by carrying out the calculation by taking into 
account the expressions for the derivatives of the unit vectors 
ei given in 3.3., or by applying the formula which gives the 
divergence of a vector field in curvilinear coordinates. The 
equations which correspond to (7) with ranks of 0 and 1 respective- 
ly are given in (17a) and (18a) . 

4.1.2. Momentum equation /35 

Here A is a tensor of second order 

A - pi + P VV + ct . 

and in order to separate B and |_ , it is first useful to expli- 

citly write down the components of the dissipative terms. The 
formulas (3-26) show that using the internal variables we have 


* 
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With 



dw 

et u = \ + e A 

dC 


,( h *Ul).» + ( h l U 2),2 


h J h 2 


+ W 



+ 2 eJT 


/ ®i.i u 2 h 1>2 K i w \ 

( — + — ) * 
\ V, h h ' 


K, w 
h. 


_ dW _ ( + < h I U 2>,2 _ / K X K 2 

et 22 = A + eJ. < + w ( - 3 — + 3 — 

dC ( hjh 2 h 2 


+ 2 ep 




2,2 U l h 2,I 

+ + 


K 2 W \ 
T~' ’ 


et. 


_ _ dW _ UMlJ.l + ( h l U 2),2 

(A+ 2 h) + e A < + w 

dC ( hjh 2 



( 10 ) 


eT. 


- €t 


13 


du i _ / w ,i 

M + « |i 

dC h 2 




this, we then have /^6 
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•f-e,.L = (c 1 - 


r i) e, + (c. 


e 2 + (C, - 


V,) e 3 


( 11 ) 


The and are primitives of the expressions, which are 
easy to write down explicitly in general, but here we will only 
give them for order zero, using a property which will be esta- 
blished in paragraph 4.2, that is that u^°^ and u 2 ^ 0 ^ are 
independent of £ and equal to and the values of u^ and u 2 

o.n the shock wave for a perfect fluid. 

Then we find: 


* 
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_ (o) (H l U 1 «J (o, ) <1 + (H^U, ? M ) it 


( K » K *\ 

+ ( 2 — + ) pw U, + 

' XJ XT ’ 


c i e 


h,h 2 


H, H 2 


p (o) 


/U.H,,, - U 2 H 2>1 , py 

u . ( ) + • 

V H.H, ' H, 


P (o) 


h 2 h 2 




_ Co) 
C, 


(H 2 U 2 U 2 p M ) 2 + (HjU* p M ) tZ Kl K 2 . 

+ / + 2 ) pw'^U, 

h 2 h 2 x h 4 h 2 7 


+ p (o, u 2 


u 2 h 2>1 - u 2 h m p™ 

+ 


h 2 h 2 


H. 


( 12 ) 


(o) + (H^pw 10 ^ /Kj 

C, + 


;<o>\ 


\ 


h 2 h 2 

_ (0) . k 2 u; 

x H 2 h 2 


/ K, K * \ 

(r + ;r) 


p w 2 


l(o) 


H, H 2 


)- 
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(13) 


r 
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We find that in (12) , considering the properties of the inter- 
nal structure of order zero which will be discussed in paragraph 
4.2, P"w (o) is independent of £ and equal to the value taken 

on by p w on the shock wave in a perfect fluid. 

The equations corresponding to (7) with respective ranks of 
0 and 1 are given by (17b, c,d) and (18b, c,d) . 


4.1.3. Energy equation 


In this case, A is a vector 

v* 

2 


/ V*\ 

( h + — j pV + e(t.V + q ) . 


(14) 


q is given by an expression in terms of internal variables as 
follows 

_ / e i _ • e. 


__ dT 

€ q = k — e 3 
dC 


- / 1 - - \ 

+ ek (=- T .. + =- T .. 

\ h, h 2 


We now have: 


( V 2 \ dw 

h + — j pw - (‘X+2P) w . - 


/ du i _ \ _ 6T 

k — + u 2 — I - it— . 

V dc dC > dC 


dZ 


6T 


. __ + ( h l U 2>,2 _/ K l K *\ 

e 3 * L = d 4 = - A.w — \T~ Y - ' " 


h 2 h 2 


1 "2 
w, K,u, 


/ K i u i \ / w - 2 ^ 2 2 \ 

- p u, — r— - p u 2 I — 

\ h. h. ' ^ h. h 2 ' 


r -e ,.L= c. 


(15) 
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and C^, are the primitives with respect to z, of the expression 
which is easily written down explicitly in general. However, we 
will only give it here up to order zero, always taking into 
accou 
of C : 


account that u|°^ = and u 2 °^ = U 2 and U 2 are independent 


/ 


ci o) 


(u,H 2 P (h + ^)' 0) ) #1 + (u 2 H lP (h + ^) <0) ), 


Hi h 2 


U, h 2 j 


v z 

P w (h + — ) 


(o) 


_ w ♦ (wiT)- 


(16) 


H, H 2 


(-l; + H;)( k 5T + (a. + 2 P ) w ) 


The equations corresponding to (7) with respective ranks 0 
and 1 are given by (17e) and (18e) . 

4.2. STUDY OF THE INTERNAL STRUCTURES OF ORDER 0 AND 1 


According to the preceding results, the internal structure 
of order 0 is defined by the nonlinear system. 


a) 

'=) 

=) (17) 


pv/ (o) = 

M ot 

> (0) 

dll. 

MoU, (o) 

- • 


dC 

— ( o) 

dii. 

Mou'”’ 

- jl <0 ’ 

dC 


_ p (°> 

i 


-<o) 


-M 


■ (o) 


\ M 0 (h+^i) - fTUj 0) 


p(o) 

2 

r 


6w (o) 

— 

= p (o) 

3 * 


dC 



— ( o) 

du, 

“a 

— (c) *“« 

, . dw <o) 6T ( 


-pu 2 - 

(A+ 2 M) w (o) -k 

d C 

d C 

dC d C 


o) 


d) 


_.te 
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in which m 6 , p°, p 2 °, Pj 0 , E 0 designates the constants. According 
to the lemma of the secohd part, the unknown functions 
p <0> , "w w , u <o) have to strive towards the finite limits when 

5 increases indefinitely in. absolute value. This condition, 
considering (17b) and (17c) , implies that u|°^ and u^ 0 ^ remain 
constant: the tangential components remain invariant and equal 

to U 1 and U 2 in the internal structure of order 0 . The problem 
is then reduced to solving a system of two equations, eseentially 
(17b) and (17c), with two unknowns, for example, p°, T°. The 
other thermodynamic quantities p°, JT, X°, ... are derived from them 
using the state equation and the law of variation of dissipative 
coefficients as a function of temperature. w° is derived from 
(17a) . This problem is a classical one [2]. One finds that for 
a system of values of the constants (M, P, E) which insure the 
possibility of a shock wave, there is one solution and only one 
solution (naturally defined up to a translation in ?) and that 
the various values of the dynamic and thermodynamic quantities 
for £ = + ro and C = -°° will satisfy the zero order shock conditions. 
This is as it has to be according to the theorem of the second part, 
section 2.3. One also can show that the functionsp ^ ( £) , 

T ( 5 ) obtained in this way have the properties of the first term 
of the internal expansion of a function of class NS (lemma 4 — 
formula 2 - 11 ) . 

We will now propose to study the internal structure of rank 1. 
The equations which govern this are written in the following form 
according to the preceding results: 



P W (1) + ci o) = 0 , 


du 


(« 


(i) 


P wu i . H 


FT(o) (o) Co) _ 


a C 


<5u 


+ c; 01 + d, ' - v, 


( 1 ) 


(1) " 2 , — Co) , — (o) — (o) 

P wUj - p + C 2 + D 2 - v 2 

a C 


= o , 


= 0 


a) 

b) 

c) 


(i) 


dw 


U) 


p + pw 2 - ( X + 2 p) 


c>C 


— (o) — (o) — -Co) 

+ C, + D, - V, =0, 


■ (l) 


p w 


K) 


■d) 


dw'*' “ST 
- ( X + 2 p) w — -k — 

SC sc 


a) 


-p-<0> 


( 1 ) 


( 1 ) 


( OU, OU- \ 

Uj + Uj, ) + C‘ 4 ° ) + D^ o) 

dr dr/ 


d) 


e) 


V 4 (o) = 0, 


( 18 ) 


These are formulas in which the c|°^ and the Vj°^ are primitives 
of the functions c|°^ and vf°^ / whose explicit expressions were 
written down in (9) , (12) , (13) , (16) . Because of the proper- 
ties which result from solving (17) and which were discussed 
above, these are functions of C. For large values of C, they 
have the behavior of polynomials with a degree which is at least 
equal to one. In contrast to (17), this system (18) is linear. 


because a quantity like a b c 
-(1) r( 1) -(1) 


( 1 ) 


is linear with respect to 


abE U) = ir> (o) c (1) + c"a (o) T^ 1 ’ + -5c (o) a <l) 


In order to simplify the notation, we will write down the 
system in the form 

' fl («)+(?>« o , (19) 

where ft designates a linear operator which operates on the flow 
l defined by w (1) . p (,) . p (1) . . . and (7) which is a vector with the 
components 

y. = C- o) + D ; (o) - V; (o) , (D 0 (o) = 0 . V 0 (o) = 0) 
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4.2.1. The associated homogeneous linear system 


/40 


Now we must study the homogeneous linear system associated 
with (18) . According to lemma (4) of the second part, 
w (1) , p (1) . p ( i) .... is limited by polynomials of the first degree 
at the most in C for large values of n which are positive or 
negative. Consequently, only the solutions of the homogeneous 
system which satisfy this property have to be retained. 


Let w, p, u ... be a solution of a ( £ ) = °* 
Of necessity 


— ( o) p + w p £ 


bw 


Uj = 0 , u 2 = 0 . 

^ dw (o> 


p + M 0 w - (A.+ 2 (i) - (X + 2 — — " 0 • 

bn &C 


■<o> 


dw 


b w 


(o) 


M o (h + w (0> w) - ( A. + 2 p ) w ( X + 2 p ) w 

b c on 


dw 


(o) 


dT 


(o) 


- ( X + 2 (i ) w -k 

bn at 


6T 

i>n 


- k' 0 * — = o . 


( 20 ) 


and consequently the system to be studied is of the form 

~ = a yf + P T , yL= r w + St , ( 21 ) 

where a, 8, y, 6 are functions of n defined and limited for all 
values of £. 


Any solution of (21) is defined for any value of C • We assume 
here, which can be verified rigorously, that its behavior for 
values of n with a large absolute value is obtained by consider- 
ing the behavior of the solutions of the linear homogeneous systems 
with constant coefficients formed from (21) , by replacing a, 8, 

Y, 6 by their values taken on respectively for £ =+°° and £ =-°°. 

If one considers (20) , it is clear that the systems with constant 
coefficients are precisely those which one would obtain when one 
would study the behavior of the integrals of (17) for C = +°° and 
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C =-°°. In the classical theory of the internal structure of order 
0, one establishes that for z, =-°° all of the solutions are expon- 
entially decreasing. For T, =+°°, the solutions in general are 
exponentially increasing with the exception of the single family 
2 of solutions which are exponentially decreasing, up to a multi- 
plication factor. We therefore can state the following: 

Lemma 1. The system (20) has a solution (?), defined up to 
a multiplication constant which decreases exponentially at infinity. 

Any other solution is exponentially increasing. 

4.2.2. Solution of the nonhomogeneous system (19) 

Considering the previous result, it is sufficient here to 
study a particular solution. Each of the components 7 t of the 
vector (7) can be written in the form 

^(^(fri) +p ‘ (S) ) + t( p+ (^) - P-(5i) ) tgh C +?. . ( 22) 

P + ( 7 ; ) and are polynomials of degree 1 and 7 is an expon- /41 

entially decreasing function of K. Let a + and a" be the differ- 
ential operators with constant coefficients obtained by setting 
Z, =+°° and £=-°° in the coefficients of & . One then insures 

that the systems 

a + ( z) = p + (y) , a" (?) «= p- (y) , ( 23 ) 

have solutions whose components w, p, p...are polynomials of degree 

1 in C. We will call these solutions P + ( Z ) and p- ( z ) which are 

* 

obtained by substitution and identification . If we now set 


+p_ i »)*t( 

P + (?i) - P' Ui) 

) tgh C + Z; 

(24) 

the system of the Z-, satisfies 

a (?) + 

the symbolic 

5 = 0 

equation 

(25) 


and it is easy to see tlfat S is a function of C which is exponen- 
tially decreasing at infinity. Here we will assume without a 
— - 

One is led to linear systems whose determinants are not zero 
when the velocity is not sonic, a condition which is obviously 
excluded. 
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proof the following lemma: 


LEMMA 2. The system (25)»in which s is a given function which 
decreases exponentially at infinity, itself has at least one 
exponentially decreasing solution at infinity. 

We can now collect the preceding results and state the fol- 
lowing : 

THEOREM. The differential system which defines the internal struc- 
ture of rank 1 has a family of solutions which depend on an arbi- 
trary constant. These are solutions which behave like polynomials 
of degree 1 in ? for large values of positive and negative z; . 

We should note that the preceding statement assumes that(T) 
is well determined. In fact, we have seen that each component 
of (^) is only defined up to an additive constant. Therefore, the 
most general solution of the internal structure of rank 1 corres- 
ponding to a given zero order structure depends on 6 constants . 

4.3. RELATIONSHIP WITH THE EXTERNAL FLOW — SHOCK CONDITIONS 

In order to provide coherence of the results, it is necessary 
that the solution found in this way has the properties of the 
terms of rank 1 of the internal structure of a system of NS func- 
tions. For this purpose, it is necessary to verify that the coeffi- 
cients of c in each of the components w :i) , ?r tl) ,p (l1 , d e p + (ZT) 
and P~ ( ; ) are identical to the derivatives ^L (o) ^£ (o) J>p (o) 

d z ’ dz ‘ dz ’ ‘ * * 

taken for z = + 0 and z = - 0 of the corres- 
ponding quantities w <o) . p<°> , p<°) , ... in perfect fluid flow. This 
proof results immediately because the equations of perfect fluids 
are the following in symbolic form 

dA M _ / ei aA < 0) e ; . aA ( o) \ (26) 

e 3 ' dz ' V h x * dx t h j dx 2 J 
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and 
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In succession A <0 ' is equal to (pV) (o) , (p) <0> 1 + (P V V) 

/ V 2 \/ \(o) 

v 11 + ~2 ) P V / ; • But the values of the second terms of (26) for 

z = +0 are precisely the values of -c^ 0 ’ for C= +°°. One can easily 
see that the first terms of (26) are identical with the product of 
the coefficients w (,) , p (1) , p (1) , ... in<a±,'j) and the coefficients of 

£ in the polynomial part of the same w (1) , p (1) , p (,) for z = + 0 . 

This result follows from the fact that the coefficients of the com- 
ponents of Pt ( s ) in (26) are precisely the c j (o) for C ± °°. Thus, 
for the mass conservation equation (26) leads to the following 




+ (C 0 ) 


= 0 




(26*) 


and (23) leads to 


/ (o) j|w_ 

V P dC 


-( 1 ) 


+ w ( °> 






+ (C 0 (0) ) 


C=t. 


(23’) 


naturally we find that (P (o) ) z= ± 0 = (P (o> )^ = ±« > as well as the relationship 
similar to w (o) . 


If we are now interested in the constant terms in each of the 
components w (1) , ... of p+ ( j ) and P ~(e) , we are led to finding the 
shodk conditions using a different method than the one used in the 
third part, that is, one which agrees with the hypothesis used in 
this third part. The various functions which characterize the 
characteristics of the flow are NS functions and this is not contra- 
dictory, and this is true up to at least order of 1. We have veri- 
fied the fact that the same thing holds true up to order n for arbi- 
trary n. Here we will be content with giving this proof for the 
fourth equation (18) only for the terms of rank 1. The procedure is 
simple: the first term of the equation is uniformly limited in C 

if one isolates a term of the form £ (a + b thO and the form of a 
and b is exactly the one which we have discussed. Once we have 
subtracted this, let us then subtract the values of the two first 
terms for £ = + 00 and ? = -°°.for each term. It is clear that 
p + pw J,1) in this way gives a result which is none other than 
[p+pw 2<1) ] . As for the term (X + 2 p) ■— 1 , the conditions are not as 

clear but if we write 
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-Jf-U) (o) XT (l) U) dw (o) 

(X+2H)|^ '(X+2|i) Yl + {X+2M) Tl 


-t-M 

and note that is zero at infinity and also what we have 

said about the terms in £ in P±(£), we find that = fpL {o) '\ 

5^(0 Uc J c . /z-*o 

We see that the term (X+2 jj)-^ gives a result 

equal to £(X + 2p)|^J > using the rule defined at the beginning 

/ \ 

of this paragraph. The term D- obviously gives [D_] , and we 

j f \ -3 #0* 

only have to examine what the result from and is. 

quantities ^ and are primitives of c^°^ and v^°^ . 

is clear that the result has to be c and v . If we refer to 
(12) and (13) taking into account that p w^°^ is independent of C 

'k 

and, therefore, gives zero contribution due to the operation we 


The 

It 


find that 


- C 3 + Vj = 



/ K » 


/ dw \ , t . 



2 p - (P w 2 ) 




V bt 1 


/K x u* 

+ p ( — + — ) 

' H , H, > 


K 2 U i' 


dw . 

= | P - A. 

6C , 


Kj_ + K«_ 


H. 


H, 


\ 


k 2 u |v 

) + P 

+ 

/ 

V H. 

H, / 


/4 3 


when we use (17d) . 

If we collect these various results, we find the third condi- 
tion (3.30) . 

In summary, the solutions discussed when explaining theorem 
4.2 satisfy all of the desired matching and coherence conditions 
with results obtained up to the order considered. 


4.4. CONCLUSION AND FINAL REMARKS 


By specifying a perfect fluid flow, which allows the presence 

of shock wave (E), one can determine the internal structure of 

* r 

order zero . From this result, it is possible to then without 

* 

difficulty write down the shock conditions of rank 1. If one 
knows the external flow of rank 1, the quantities which define 
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this flow on (I) will verify the rank 1 shock conditions. It is 
then possible to find a solution of the internal structure of rank 
1 which is compatible with the given external flow. But this solu- 
tion is not unique. It depends on an arbitrary constant. The 
difference between the two solutions decreases exponentially in t , . 

It is then clear that the shock conditions of rank 2 depend on this 
constant, even though the explicit calculations were not given here. 

This allows one to assume that this indetermination of the 
internal structure of rank 1 is undoubtedly limited to an over- 
determination of the external flow of rank 2. But this is a sug- 
gestion which seems difficult to prove. The verification will then 
involve the properties of a global problem and goes beyond the 
scope of this article in which we only consider the local problem. 

This does not affect the essential point of the results 
obtained nor their coherence. The shock conditions and the equa- 
tions for the shock structure for any rank n, at least in principle, 

have been developed. The complete results were given explicitly 
* * 

for rank 1 

NOTES ADDED UPON CORRECTING THE PROOFS 

* The internal structure of order 0 is in fact determined up 

to a translation in £ (see system 4.17) . The effect of this trans- 

★ 

lation obviously affects the evaluation of quantities such as p , 

★ 

p ,.... Therefore, there is an indeterminancy (depending on the para- 
meter) in writing down the order 1 shock conditions. This indeter- 
minancy comes from the fact that we have specified the position 
of the shock at z = 0 . In the global problem, the position of the 
shock is not arbitrary. , When writing down the shock conditions, it 
is usual to find an indeterminancy which corresponds to the over- 
abundance of conditions relative to the boundary problems which 
govern the external development. 
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** Since writing this manuscript in December 1960, the results 
of this work have been expanded to the case of non-steady flow 
and shock conditions were given explicitly using intrinsic nota- 
tions (C.R.A. Sci. t.252, p. 1101-1103, February 20, 1961). Since 
this time, we have become aware of the memoir by R. R. Chow and 
L. Ting: Higher order theory of curved shock, JAS 1961, p. 428, 

which discussed the same problem using methods which are very 
close in their technical details. 


Manuscript shipped December 1, 1960 
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